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Abstract 

When a continuous symmetry is spontaneously broken in nonrelativistic theories, there appear 
Nambu-Goldstone (NG) modes, the dispersion relations of which are either linear (type I) or 
quadratic (type II). We give a general framework to interpolate between relativistic and nonrela¬ 
tivistic NG modes, revealing a nature of type-I and -II NG modes in nonrelativistic theories. The 
interpolating Lagrangians have the nonlinear Lorentz invariance which reduces to the Galilei or 
Schrodinger invariance in the nonrelativistic limit. We find that type-I and type-II NG modes in 
the interpolating region are accompanied with a Higgs mode and a chiral NG partner, respectively, 
both of which are gapful. In the ultrarelativistic limit, a set of a type-I NG mode and its Higgs 
partner remains, while a set of a type-II NG mode and its gapful NG partner turns to a set of 
two type-I NG modes. In the nonrelativistic limit, the both types of accompanied gapful modes 
become infinitely massive, disappearing from the spectrum. The examples contain a phonon in 
Bose-Einstein condensates or helium superfluids, a phonon and magnon in spinor Bose-Einstein 
condensates, a magnon in ferromagnets, and a kelvon and dilaton-magnon localized around a 
Skyrmion line in ferromagnets. 


1 


I. INTRODUCTION 


When a continuous symmetry is spontaneously broken in nonrelativistic theories, there 
appear Nambu-Goldstone (NG) modes, the dispersion relations of which are either linear 
(type I) or quadratic (type II). The numbers of type-I and -II NG modes satisfy the Nielsen- 
Chadha inequality jlj|. After the crucial observation Q, the numbers of type-I and -II NG 
modes were summarized as the Watanabe-Brauner relation between those numbers and the 
rank of a matrix consisting of the commutation relations of broken generators evaluated 


in the ground state [sj]. The relation classifies types A and B instead of types I and II. 

"•‘TT *'■" 


method jb] and later by the Bogoliubov theory 

In the presence of topological solitons or defects, there appear NG modes localized around 
them such as translational zero modes. When NG modes are normalizable such as a domain 
wall |8] and Skyrmion line jb, 10] in ferromagnets, localized type-B NG modes have usual 
quadratic dispersion relations and are of type II. On the other hand, when NG modes are non- 
normalizable such as a domain wall in two-component Bose-Einstein condensates (BECs) 
and a vortex in scalar BECs or 1 1 Ie superfluid, type-B NG modes have usual quadratic 


dispersion relations and are of type II when the 
with the wavelength of NG modes (see, e.g . Refs. 


relations in infinite system sizes 


a 


12 


ransverse system size is small compared 
111]), but they have noninteger dispersion 
13]. The formulas of dispersion relations interpolating 


small and large system sizes were recently obtained in Ref. 


14]. A relationship between the 


number of NG modes and the homotopy group for topological solitons was also studied 
Other developments include, for instance, space-time symmetry breaking 


15|- 


metry breaking accompanied with the Higgs mechanism ITT 191. hnite temperature and 


density 


20], topological interaction 


2l| . and quasi-NG modes 


16], gauge sym- 


22 ], 



Parameters 

Symmetry 

Type-I NG mode 

Type-II NG mode 

Ultrarelativistic 

p —>• 0 

Lorentz 

1 type-I + 1 Higgs 

2 type-I 

Relativistic 

0 < c, /j, < oo 

Lorentz 

1 type-I + 1 Higgs 

1 type-II + 1 gapped 

Nonrelativistic 

c —^ oo 

Galilei (Shrodinger) 

1 type-I 

1 type-II 


TABLE I. Interpolation of type-I and type-II NG modes between nonrelativistic and ultrarelativis- 
tic theories. 
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In general, it is usually said that only type-I NG modes are possible in Lorentz invariant 
theories. It is, however, not yet clear how NG modes are interpolated between relativistic 
and nonrelativistic theories, in particular, how type-II NG modes in nonrelativistic theories 
reduce to type-I NG modes in relativistic theories when both the theories are interpolated. In 
this paper, we clarify how NG modes are interpolated between relativistic and nonrelativistic 
theories, as summarized in Table Q] We first consider relativistic Lagrangians and introduce 
a chemical potential for particles. The resulting Lagrangians, containing both first- and 
second-order time derivative terms, interpolate relativistic and nonrelativistic Lagrangians 
in the two limits: the second time derivative vanishes in the nonrelativistic limit c —> oo with 
the speed c of the light, while the first time derivative vanishes in the relativistic limit /i —?• 0, 
in which the chemical potential is sent to zero. The latter case is often called ultrarelativistic, 
so we use this terminology because the Lorentz invariance exists in the intermediate region. 
We first point out that interpolating Lagrangians have the nonlinear Lorentz invariance, 
which reduces to the Galilei or Schrodinger invariance in the nonrelativistic limit c —>■ oo. 
The Lorentz invariance becomes manifest in the ultrarelativistic limit /i —>■ 0. We find that 
there can exist either a type-I or type-II NG mode in the intermediate Lagrangian, even in 
the presence of the Lorentz invariance. Correspondingly, the Watanabe-Brauner relation 
holds in the intermediate region even in the presence of the Lorentz invariance; a commutator 
of two generators does not vanish for a type-II NG mode. We also find that each of either the 
type-I or type-II mode is accompanied with a gapful mode. The gapful mode accompanied 
with a type-I NG mode can be identified with a Higgs mode, while that accompanied with 
a type-II NG mode can be referred as a “chiral partner” or a “gapful NG partner”. In 
the ultrarelativistic limit, a set of a type-I NG mode and its Higgs partner remains as it 
is, while a set of a type-II NG mode and gapful NG partner becomes a set of two type-I 
NG modes. On the other hand, in the nonrelativistic limit, both the Higgs partner of a 
type-I NG mode and gapful NG parter of a type-II NG mode become infinitely massive and 
disappear from the spectrum, and there remains only the type-I or type-II NG mode. This 
mechanism reveals a nature of type-I and II NG modes. We show these in typical examples 
of both bulk NG modes and NG modes localized around a topological soliton. The bulk 
examples contain a phonon in scalar BECs and a magnon in ferromagnets, while soliton 
examples contain a kelvon and dilaton-magnon localized around a Skyrmion line in isotropic 
ferromagnets. An another example is a ripplon-magnon localized around a domain wall in 
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anisotropic ferromagnets studied in Ref. [8j]. In all examples, we give two approaches: the 
effective theory and linear response theory (the Bogoliubov-de Gennes equations). 

The interpolating Lagrangians containing both first- and second-order time derivatives 
that we consider in this paper appear in various contexts of both theoretical and experimen¬ 
tal physics, and thereby our results yield suggestions of several theoretical and experimental 
works. As we denoted above, they describe relativistic held theories with the finite chemical 
potential. In the quantum mechanical framework, it is well known that these models natu¬ 
rally give the complex probabilities in the path-integral formalism as long as the chemical 


potential is finite [23, 24], and we can expect some qualitative change from zero to finite 


chemical potentials. Even in the semiclassical framework for symmetry-broken systems, our 
results show that there is a drastic change of low-energy modes both in bulk and topo¬ 
logical defects: the coupling of the type-I NG mode and Higgs mode or the coupling of 
two type-I NG modes to one type-II NG mode. The interpolation between the relativistic 
and nonrelativistic frameworks appears in various ultracold atomic systems. Recently, the 
existence of Higgs mode in strongly interacting lattice bosons close to the superfluid-Mott 


insulating transition point has been theoretically 


25 


30] and experimentally [3l| proposed 


and confirmed. In this system, the first-order time derivative term is prohibited, and the 
second-order time derivative term becomes important at the transition point because of the 
particle-hole symmetry, by which the Higgs mode can be expected. In the superfluid phase 
far from the transition point, it is well known that the first-order time derivative term is 
dominant and the Higgs mode is absent. Our results explain how the Higgs and NG modes 
are changed as the parameter changes to the transition point. We show two two examples. 


One is the charged fermionic systems close to the BEC-BCS crossover point 


32] that has 


been theoretically predicted to contain the both Erst- and second-ordered time derivative 
terms, predicting the existence of the Higgs mode. The other example is the magnetic 
model 33]. As already known, two type-I NG modes exist in antiferromagnets, while one 
type-II NG mode exists in ferromagnets. In a canted ferromagnet between the ferromagnet 
and antiferromagnet, it has been reported that there appear one type-II NG mode and one 
gapful Higgs mode, which is quite similar to what we obtain in the interpolating region 
between relativistic and nonrelativistic models even though the magnetic model itself is 
nonrelativistic. 

This paper is organized as follows. In Sec. [Til we discuss bulk NG modes. We study 
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phonons in a scalar BEC and magnons in ferromagnets as examples of type-I and -II NG 
inodes in Secs. Ill Al and III Bl respectively. In Sec. IIII1 we discuss NG inodes localized around 
topological solitons. We study kelvons and dilaton-magnons localized around a Skyrmion 
line in ferromagnets as examples of type-II NG modes. Section HVl is devoted to a summary 
and discussion. In Appendix [AJ we give a further example of a spinor BEC that contains 
both type-I and -II NG modes in the bulk. 

II. NG MODES IN THE BULK 

A. Interpolating type-I NG mode: Phonons in scalar BEC 

Let us consider the Lagrangian density for a single complex scalar field i/j, interpolating 
relativistic and nonrelativistic theories, 

C = ~ ~ IWf - |(|#' - p) 2 , (1) 

where g is the coupling constant and p is the real positive constant giving a vacuum expec¬ 
tation value. This Lagrangian density interpolates between two extreme cases; it reduces 
to the relativistic Goldstone model in the ultrarelativistic limit /i —)■ 0 and to the Gross- 
Pitaevskii model in the nonrelativistic limit c —> oo. In the generic region, the Lagrangian 
density is invariant under the Lorentz transformation 

t, = 'y( t ~-Tb x ' = l i x ~vt), da = ydt + yvd x , d,a = ^d t + r yd x , 

V C 7 C (2) 

^ = e^, S = -/m 2 {(l- 7 )f + ^}, 7 = 1 

l c 2 ) yl~ v / c 

which reduces to the Galilei or Schrodinger transformation in the nonrelativistic limit c —y 
oo. The equation of motion for reads 

+ 2 ipd t i) = -V 2 ^ + g{\i>\ 2 ~ p)i>, (3) 

c 

which has the static constant solution ij; 0 = yfp. 

1. Low-energy effective theory 

The low-energy dynamics around the static solution ^ 0 can be discussed by considering 
the low-energy effective theory. We introduce fluctuations of amplitude f(x,t ) and phase 
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9(x,t) around i/jq: 


«i = V>o{l + f(x,t)}e' el -*-‘>. 


(4) 


Inserting Eq. (jj]) into Eq. (JT|) , we obtain the effective Lagrangian density 

— = t-th - 2/x(l + 2 })» - (V/) 2 - (V9) 2 - 2 gpf + 0((/, «) 3 ). 

p c 2 

The low-energy dynamics of / and 6 derived from the Euler-Lagrange equations reads 
4 + 2 pO- V 2 / + 2 gpf = 0, 4 - 2/'./' ~ V 2 0 = 0. 


(5) 


( 6 ) 


Typical solutions of Eq. (EJ) are / = /o cos(fc ■ x — cut + 5), 6 = 9 0 sin [k ■ x — cut + §) with 
the dispersions 


= icy/ 2/i 2 c 2 + k 2 + gp + 4/i 4 c 4 + 4 p 2 c 2 (k 2 + gp) + g 2 /4 

a;± G = icy/2yU 2 c 2 + k 2 + gp — \J\ /i 4 c 4 + 4 p 2 c 2 {k 2 + gp) + g 2 4. 
In the long-wavelength limit for k —» 0, these dispersions reduce to 


= ±c-\/4/i 2 c 2 + 2gp + 0(k 2 ), cu± G = ±cfc,/ 9P —— + 0(fc 2 ) (8) 

V 2c 2 /i 2 + gp 

giving rise to one gapful (wr) and one gapless (cung) mode, identified as Higgs and NG 
modes, respectively. The amplitudes / 0 for o;± and a;^ G are obtained as 


fo± = =F 


9 Q sj2p?c 2 + gp 
\f2pc 


+ 0(k 2 ), / 0 N ± G = ± 


9 0 pck 


y/9p(^ 


2 c 2 + g 2 p 2 ) 


+ 0(i 2 


(9) 


This is the relationship between amplitudes of fluctuations for coupled dynamics of / and 
0. In the long-wavelength limit k —> 0, we obtain —> 0, indicating that the oscillation 

of / vanishes and there remains the oscillation of 9 as the pure phase mode. 

In the ultrarelativistic limit p —» 0, the dynamics of / and 9 in Eq. (JSJ) are independent 
of each other, and amplitudes /o and 9q become independent variables with dispersions 


l—> ±c\/k 2 + 2gp, o;^ G —* ±ck. 


( 10 ) 


The gapful dispersion reduces to that for /, while the gapless dispersion cn^ G reduces that 
for 9, i.e., the Higgs and NG modes get to pure amplitude and phase oscillations respectively. 
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In the nonrelativistic limit c —> oo, the Higgs mode disappears with the divergent dispersion 
u;± —y oo. The NG modes remain coupled oscillations of / and 6 with 


ng v k^k 2 + 2gp 

- Tg. -’ 

In the long-wavelength limit k —> 0, the NG 
arbitrary p. 



Opk 

Vk 2 + 2 gp 


( 11 ) 


mode is always the pure phase mode with 


2. Linear-response theory 


In the linear-response framework, the dynamics of ip can be written as ip —> ipo + u + v* 
with fluctuations u and v*. Here, rewriting u and v as u = uoe l ^ kx ~ ut+s \ v = Voe P k ' x ~ ult + s ) 
with u 0 , v o G M, and inserting ip into Eq. (J3]) leads to the Bogoliubov equation 


cu 2 /c 2 + 2p,u — k 2 — gp —gp 

—gp oj 2 /c 2 — 2pu — k 2 



+ O((u 0 ,v 0 ) 2 ) = 0 (12) 


H NG 

with the dispersion relation oj = cu±’ . The fluctuation dip = u + v* becomes 


dip oc cos(fc ■ x — + 5) + 


2p 2 c 2 ± pc^Ap 2 c 2 + 2ffp ^_ i(fc . a ._ h ,H f+a) 


9P 


± t+d) +0(k 2 ), (13) 


for Higgs modes with the dispersion and 


dip oc i sin (k ■ x - ul G t + d)± === -Hk-x-^t+s) + Q , k 2 n ( 14 ) 

y/gpQ^ + gp) 

for NG modes with the dispersion cu^ G . Because the ground state ipo = y/p has only the 
real part, the first terms of the right-hand sides in Eqs. (1T3T) and (1T41) can be regarded as the 
amplitude and phase oscillations respectively. Both the second terms of the right-hand sides 
in Eqs. (1T3|) and (1T4|) are coupled oscillations of the amplitude and phase. They vanish in the 
ultrarelativistic limit, which reveals that the Higgs and NG modes become purely amplitude 
and phase oscillations, respectively. In the nonrelativistic limit, on the other hand, the NG 
mode 


dip oc % sin (k • x - + d) ± ± G *+ 5 ) + 0(k 2 ), (15) 

V2 gp 

remains to be a coupled amplitude and phase oscillation. In the long-wavelength limit k —> 0, 
the NG mode is always the pure phase oscillation with arbitrary p. Equations. (fl3]l and dl~4]) 
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are consistent with the expansion of -0 = y/p-fl + f^± lG cos(k ■ x — u)± ,>>l(J t + 8)} exp{f6* 0 sin(fc- 


H,NG, 


X 


w^’ NG i + 5)} around 9 0 = 0, which reveals that both the low-energy effective theory 


and linear-response theory give the same Higgs and NG inodes. Similar behaviors of NG 
and Higgs modes are theoretically reported for strongly interacting lattice bosons close to 
superfluid-Mott insulating transition point 281430]. 


B. Interpolating type-II NG mode: Magnons in ferromagnets 


We consider the interpolating Lagrangian density for the continuum Heisenberg model 
or the 0(3) nonlinear sigma (CP 1 ) model 


C = 



ft 

2 

ifi{u*ii - 


ii*u ) Vw : 

2 

c 2 (l + 

\u\ 

to 1 

to I 

- 

+ 

\u\ 

to 

+ 

\u\ 

2)2 


(16) 


where u E C is the complex projective coordinate of CP 1 , defined as (j) 7 = (1, u) T /yfl + \u\ 2 
with normalized two complex scalar fields 0 = (0 i,02) T with |0!| 2 + 10 2 1 2 = 1. This La¬ 
grangian density is invariant under the following Lorentz transformation: 


t' = , x = y(x - vt), df = 7 d t + ”/vd x , d x > = d t + 7<9 X , 

u' = e lS u, d t S = — (1 — 7)/ic 2 (l + |w| 2 ), d x S =+ \u\ 2 ). 


(17) 


This reduces to the Galilei or Schrodinger transformation in the nonrelativistic limit c —> oo 
The equation of motion for u reads 

(|u| 2 + 1 )u - 2 u*u 2 


2f/i(|u| 2 + 1 )u = (\u\ 2 + l)V“it — 2 m*(V«) 2 , 


(18) 


which has the uniform and static solution w 0 = const. 

Under the Hopf map for a 3-vector of real scalar fields n = 0tcr0 with the Pauli matrices 
cr, the Lagrangian density (ITU]) describes the isotropic Heisenberg ferromagnets 


C = ~ 

c z 


\n 


/ic 2 (hi77. 2 - nih 2 ) 1 |Vn| 2 

4 1 + n 3 J 4 


+ 


(19) 


1. Low-energy effective theory 

Here, we consider the low-energy effective theory for the low-energy excitation, with fixing 
a uniform and static solution uq = 0 and its fluctuation Su = a + i/3 with a, (3 e K. In 













terms of n, uq = 0 is equivalent to 713 = 1, and a and f3 are the fluctuations of n\ and n 2 , 
respectively. The effective Lagrangian density becomes 

. 2 <o 2 

Ces = a + f + 2fi{a/3 - a(3) - (V « 2 + V/? 2 ) + 0((a, p) 3 ). ( 2 0) 

cr 

The low-energy dynamics of a and f3 becomes 

% + 2^/3 - V 2 a = 0, 4-2/id - V 2 /3. (21) 

cr cr 

As in the previous case for phonons in a scalar BEC, the dynamics of a and f3 are independent 
of each other only in the ultrarelativistic limit /i —y 0. Typical solutions are a = a 0 cos (k ■ 
x — cot + 5), f3 = (3 q sin(fc ■ x — cut + 5), with 


to. 


CO. 


NG 

± 


= ±c^y / fc^+7i 2 c 2 + l^cj = ±^2/ic 2 + + 0(k 4 ), ao ± 

= ±c(\/k 2 + jj?c 2 - ficj = ±^~ + 0(k 4 ), a^ G = ±@q±. 


iH 

0 ±> 


( 22 ) 


The second solution is a type-11 NG mode which is a magnon, and the first one is its 
chiral massive partner which we may call a “massive magnon”. In the ultrarelativistic limit 
H —> 0, Eq. ()22l) reduces to u;±’ NG = ±ck and cko± N& an d /^o± NG are independent of each 
other, which implies that the type-II NG and Higgs modes reduce to two type-I NG modes. 
In the nonrelativistic limit c —>■ 00 , on the other hand, it reduces to 


—> 00 , 


00. 


NG 


k 2 

-> ±—. 
2/i 


(23) 


While the type-II NG mode remains gapless, the chiral massive partner becomes infinitely 
massive and disappears from the spectrum. 


2. Linear-response theory 

We consider the dynamics of magnons in the linear-response theory framework: u = 
a + e i{k-x-uit+5) _j_ a _ e -i(k-x-ut+6 ) w j^} 1 a± (= Inserting this ansatz into the dynamical equa¬ 
tion ([HD , we obtain the Bogoliubov equation 

± 2/ico^J a± = k 2 a± + 0(a|), (24) 

giving the dispersion relation 00 ^ and o;^ G in Eq. (J22]) with arbitrary a± ^ 0. The gapful 
mode for a ^ e ^h k - x ^± t + 5 ) anc j yjg mode for a± G e ± *( fc ' iE_a, ± Gt+<5 ) p r() p a g a t e i n the direction 
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parallel to k for the upper sign and antiparallel to k for the lower sign, respectively, where 
their chiralities are opposite to each other. In the ultrarelativistic limit, = o;^ G = ±ck 
leads a^{e T *( fc ' a: - a, ± t+ ' 5 ) + e ±i ( fc -*-^± G i+< 5 )| — 2a^ cos(fc • x =|= ckt + 6) with = a± G and 
a ^ e : f l ( kx -^± t + s ) — e ±*( fc - a; -^± G i+< 5 )| — =p2ia^sin(fc • x =|= ckt + S ) with = — a± G , which 
give purely real and imaginary modes. These results completely agree with those obtained 
in the low-energy effective theory. 


III. NG MODES LOCALIZED AROUND SOLITONS 


A. Kelvon and dilaton-magnon of a Skyrmion 


We start from the CP 1 Lagrangian density £ in Eq. (fl6|) . Instead of the uniform u, we 
consider a straight Skyrmion-line solution 34] extending along the z axis, 


u B (x,y,z) = 


re 


i(e+e) 


R s + R 


, r = \/{x — X ) 2 + (y — Y) 2 , 9 = tan 


-i y-Y 
x — X' 


(25) 


as the static solution. Here, R s G M + is the characteristic size of the Skyrmion line, and X, 
Y G M, 0 < 6 < 27t, and R G R are the translational, phase, and dilatation moduli of the 
Skyrmion line, respectively. 


1. Low-energy effective theory 


To discuss the low-energy dynamics of the Skyrmion line, we use the moduli approxi¬ 


mation 


351 ]; we introduce the z and t dependences of these four moduli and integrate the 


Lagrangian density in the rcy-plane with radius L : 

X 2 + Y 2 


pL py/L 2 —x 2 

I dx I dy jC = —27r + n 

' —L J-VL 2 -x 2 


2ti log 


R 2 + R 2 e 2 


c 

-)2q2 


(■Xl + U 2 ) + MXY - XY )} 


- (It: + lt:0:i + 2^9 + 2 R,R0) ^ + 0((X. y. e, Rf) 


(26) 


Rsj l cf 

The first 2n term in the right-hand side is the tension (the energy per unit length) of the 
Skyrmion line. The low-energy dynamics of X, Y, 9 , and R derived from the Euler-Lagrange 
equation becomes 

X = c 2 X zz - 2/ic 2 Y, Y = c 2 Y zz + 2 nc 2 X, (27a) 

R s 9 = c 2 R s 9 zz — 2 /j,c 2 R, R = c 2 R zz + 2 yc 2 R s 9. (27b) 
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Equation (I27a,p has the same form as Eq. (12T| with rewriting V —> d z , aq —> X, and 
«2 —> Y. Typical solutions X = X 0 cos {kz — cat + 5) and Y = Losing;? — cut + S) are 
therefore the same as those in Eq. ([22])• As long as /i ^ 0, the two translational moduli 
X and Y couple to each other forming gapless and gapful helical modes with X 0 = Y 0 
and X 0 = —Y 0 propagating along the z direction helically and antihelically. The former is 
nothing but a helical Kelvin wave or a helical kelvon if quantized as a particle, while the 
latter may be called a “massive helical kelvon.” The moduli fields 9 and R have the solution 
9 = 9 0 cos (kz — cut + 5), R = Ro sin(/cz — out + 5) with the dispersion shown in Eq. (1221) . 
The phase and dilatation moduli 9 and R couple to each other forming gapless and gapful 
modes with 9 0 = Rq/R s and 9 0 = —R 0 /R s propagating along the z direction helically and 


antihelically. We called the former a “dilaton-magnon” 


10], and the latter may be called a 


“massive dilaton-magnon.” In the ultrarelativistic limit 0, the four modes for A", Y, 9, 
and R propagate independently of each other with the linear dispersion ck as wavy kelvons 
for X and Y, U( 1) rnagnon for 9 , and dilaton for R. 

Here, we note that the dilatation symmetry is not the symmetry of the Lagrangian density 
(HBD but the symmetry of the stationary state of the dynamical equation for the Lagrangian, 
and the dilaton is not the NG mode but a so-called quasi-NG (QNG) mode [22], while 
wavy kelvons and phonons are NG modes. The dilaton-magnon is also regarded as coupled 
NG-QNG mode, while a helical kelvon is a coupled NG mode. In the nonrelativistic limit 
c —y oo, the massive helical kelvon and massive dilaton-magnon disappear because of the 


divergent c 
In Ref. 


ispersion relation. 

8], NG modes localized around a domain wall in ferromagnets with one easy axis 


were studied. The model is a nonrelativistic version of the massive CP 1 model often studied 


in the supersymmetric context 


The domain wall breaks the translational symmetry 


transverse to the wall as well as the internal 1/(1) symmetry. There appear an associated 
ripple mode and 17(1) NG modes, coupled to each other. The interpolation between ultra¬ 
relativistic and nonrelativistic theories is parallel to the case of a Skyrmion line. 


2. Linear-response theory 

We consider the dynamics of the Kelvin wave and dilaton-magnon in the linear-response 
theory framework: u = u s (R = 9 = X = Y = 0) + a+e l ( k - x ~ ut + s ) + a _ e -i(k-x-wt+5) _ i nser ting 
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this ansatz into the dynamical equation (fT8j) . we obtain the Bogoliubov-de Gennes equation: 


u 


± 2yu )a± = <(k — V;) + 


4 (rd r ± idg) 
r 2 + R 2 


a± + 0(u±), 


(28) 


where, V r = (d x ,d y ) denotes the derivative in the xy plane. By expanding a± as a± = 
a ±,i e%ie i we obtain 


or 


— ± 2/rcu a±j = l (k - d r - d r /r + l /r 2 ) + 


4(r<9 r =F /) 
r 2 + R 2 


a ±,i + 0(a±). (29) 


There are two characteristic solutions: a-po oc 1 with / = 0 and a±,\ oc r/R s with l = 1 with 
the dispersion relation shown in Eq. ()22il . As long as /i ^ 0, there are a gapless NG mode 
with o;^ G and gapful Higgs mode with cu±, and the solution becomes 


?/ NG - 
u ±,o ~ 


re 


w 


R s 


X oe ±i ( kz ~ UJ ± Gt + s '> ; 


<0 = 


re 


w 


R* 




for / = 0 and 


NG _ re ie i6 0 re ±i{ - kz - u ± Gt+5) 
“ "i“ 


«? = 


re ie id 0 re Ti{ - kz -^ t+s) 


+ 


(30) 


(31) 


R s R s R s R s 

for / = 1. w(j? G , w± 0 , and u± i are equivalent to the solution (J25]) with moduli X , Y, 
9, and R for helical kelvon, massive helical kelvon, dilaton-magnon, and massive dilaton- 
magnon, respectively, in the first order of Xq and 9q. In the ultrarelativistic limit y —y 0, 
both a;± G and thave linear dispersion relations ±ck, giving wavy kelvons as the linear 
combination of u^ G and u± 0 and the 1/(1) magnon and dilaton as the linear combination 
of m± g and u® v 

We shortly note other solutions having the same dispersions o;^ G and With l = 0 
and l = 1, solutions a± ]0 and a± tl have their anomalous pairs a± )0 oc ( r/R s ) A + 4 {r/R s ) 2 + 
41og(r/i? s ) and a-pi oc ( r/R s ) 3 — ( R s /r ) + 8 (r/R s ) log (r/R s ), which are the modes changing 
the Skyrmion charge of the total volume in which we are not interested. For higher l > 2, 
there are also solutions a±,i >2 oc (r/ R s ) 1 . They do not change the state around the Skyrmion 
at the center but change the bulk state far from the Skyrmion, giving bulk magnons with 


u;± G and a bulk massive magnon with propagating along arbitrary directions. For lower 
l < 0, solutions GLtjico oc (Rs/ r )~ l give the Skyrmion-splitting modes from 1 Skyrmion at 
the center with the charge +1 to 1 Skyrmion with the charge —l at the center and (Z + 1) 
Skyrmions with the charge +1 around the center. As well as the dilaton, these Skyrmion- 
splitting modes do not come from the symmetry of the Lagrangian density, and can be 
regarded as QNG modes and their massive partners. 
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IV. SUMMARY AND DISCUSSION 


In summary, we have revealed how relativistic and nonrelativistic NG modes are inter¬ 
polated. We have found that type-I and type-II NG modes in the interpolating Lagrangians 
with the Lorentz invariance are accompanied with a gapful Higgs mode and gapful chiral NG 
partner, respectively. In the ultrarelativistic limit, the type-I NG and Higgs partner remain, 
and the type-II NG mode and gapful NG partner become two type-I NG modes. In the 
nonrelativistic limit, the accompanied gapful modes become infinitely massive, disappear¬ 
ing from the spectrum. In the whole region, the commutation relation holds consistently, 
showing that the Lorentz invariance does not forbid type-II NG modes. 

While we have studied a kelvon localized around a Skyrmion line in ferromagnets, we 
have not studied a kelvon localized around a vortex line in scalar BECs. In the latter 
case, the dispersion relation of the kelvon depends on the transverse system size. When 
the size is finite, the dispersion relation is quadratic, but it is not an integer anymore for 
infinite system size. Recently, the interpolating formula of the dis per sion relation of the 


kelvon for an arbitrary transverse system size was obtained in Ref. 


14]. Extending to the 


case of a kelvon of a vortex in an arbitrary system size for interpolating relativistic and 
nonrelativistic systems is an interesting future work that would be important for the Mott 
insulator transition point in BECs in the optical lattice. 

One of the related topics is the relaxation dynamics of topological defects generated 
through the Kibble-Zurek mechanism after the temperature quench. It has been predicted 


37] that the relaxation dynamics is universal and dependent only on a small number of 


factors such as conserved quantities, external currents, viscosity, and off-criticality. For 
one of the future problems in this topic, we can consider the dependence of the relaxation 
dynamics on the dispersion relation of NG modes excited along the topological defects in 
the interpolation between relativistic and nonrelativistic regions, which is expected to be an 
essential problem in the phase-ordering dynamics of U( 1) bosons with the finite chemical 
potential [ssj]. 

Finally, let us mention quantum corrections of NG modes in lower dimensions. In nonrel¬ 
ativistic limit, type-II NG modes remain gapless under nonperturbative quantum corrections 


391 ]. as opposed to type-I NG modes that become gapful to be consistent with the Coleman- 


Mermin-Wargner theorem. It is interesting to see whether quantum corrections give gaps 
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to type-II NG modes in the intermediate region. 
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Appendix A: Spinor BEC 

Here, we discuss a ferromagnetic F — 1 spinor BEC which contains both type-I and -II 
NG modes simultaneously. We see that it contains further a pair of gapful modes. The 
interpolating Lagrangian is 

£ = \\d t i>\ 2 + zp('0 t <9t'0 - dtipty) - |V0| 2 - ^(|-0| 2 - ——^-X (Al) 

t 2 ( g 0 J 2 

where i/j = ("0i, V’O) '4 > -i) T is the three-component (spinor-1) complex scalar fields and F 
is the triplet of the 3 by 3 50(3) generators (spin-1 spin matrices). Besides the Lorentz 
transformation given in Eq. (J2]) for all components "0±i and 0 O > this Lagrangian is invariant 
under the shift of the overall phase 0 —t 0e* e and the 50(3) spin rotation 0 —y ^>e~ lF ' s . 
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When the two coupling constants go and g\ satisfy go > g\ > 0, there exists a stable and 
static solution i/j g = (yfp, 0, 0) 7 for F = = (0, 0, p ) as the ground state. 

As well as the previous examples, we consider the following low-energy excited state: 

= Vp(( 1 + /i)e* 01 , «o + ifo, «-i + i/3- i) T , (A2) 


where /i and 9 i are fluctuations of the amplitude and phase of the first component of i/j 
and a m and /3 m are the real and imaginary parts of fluctuations of the mth component 
(m = 0, —1) of t/j. Inserting Eq. (IA2|) into Eq. (I A1 [) . we get the effective Lagrangian 


£ _ fi + 9\ + + Po + ®-i + /3-i 

p c 2 

- 2/i{(l + 2fi)9i + «o/3o — do/3o + cc_i/3_i — b-i/3_i} 

- (|V/i| 2 + |V^| 2 + |Vcc 0 | 2 + |V/3 0 | 2 + |Vcv_!| 2 + iV/Ail 2 ) 

- 2(p 0 - 9i)pfi ~ 2g l p(a 2 _ 1 + /3^) + 0((/i, 0i, a 0 , A>, a_i, /3_i) 3 ). 


(A3) 


The low-energy dynamics becomes 

^ + 2/A - V 2 /i + 2(p 0 - gi)ph = 0, j 2 - 2/A - W 2 9 1 = 0, (A4a) 

^ + 2p$ 0 ~ V 2 a 0 = 0, ^ - 2pd 0 - V 2 /3 0 = 0 , (A4b) 

c c 

—^—|- 2p/3_i — V 2 q:_i + 2gipa\ = 0, ——— 2pa—\ — V 2 /3_i + 2gipf3± = 0. (A4c) 

Equation (lA4al) has the same form as Eq. (J 6 J) with rewriting / — > fi, 9 —> 9i, and g —* go—gi- 
Therefore, as long as /i ^ 0, dynamics of fi and 9\ are coupled as fi = fio cos (k ■ x — cut + 5) 
and 9i = 9 w sm(k ■ x — cot + 5) with the dispersions shown in Eqs. (J7J) and (J9]): one gapful 
Higgs mode with n/± and one type-I gapless NG mode with cu^ G . In the relativistic limit, 
two dynamics of fi and 9\ are independent of each other giving a Higgs mode for f\ and a 
type-I NG mode for 9\. In the nonrelativistic limit, the Higgs mode vanishes with a diverging 
spectrum. 

Equation (IA4bjl has the same form as Eq. (ITT)) , and the solutions a 0 = a 0 o cos (k-x—oot+5) 
and /?o = /Soo sin(fc ■ x — cot + 6) exactly behave as a and f3\ one Higgs mode with to = 
and (3 oo = T«oo and one type-II NG mode with uo = and /?oo = ±aoo having opposite 
chiralities as long as /i 7 ^ 0. In the ultrarelativistic limit, the Higgs and type-II NG modes 
are degenerated, giving rise to two type-I NG modes. In the nonrelativistic limit, the Higgs 
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mode vanishes with a diverging spectrum. The modes a 0 and (3q can be considered as 
fluctuations of the spin rotation around F = (0, 0,p). Fluctuations for F x and F y can be 
written as e~ lFyS/ ip g = ^(1, s/\/2, 0) T + 0(s 2 ) and e lFxS iJ > g = y/p(l, is/y/2, 0) T + 0(s 2 ). 
The real and imaginary parts a o and (3q therefore correspond to fluctuations of F x and F y , 
respectively, which is consistent with a and (3 in Eq. (1201) as fluctuations of n\ and n<i- 
For Eq. (1A4cl) . typical solutions are a_i = a_i 0 cos(fc • x — cot + 5) and /3_i = (3- i0 sin (A; ■ 
x — cot + S) as long as p ^ 0 with the dispersion 


u. 


G1 

± 


LU. 


G2 


= ±c(^\/k 2 + p 2 c 2 + 2g 1 p + pcj 

/ p 2 

= ±c v 7 h 2 c 2 + 2giP + W 0 „ 

V 2y'iic 2 + 2gip 

= ±C ^ ^/k 2 ^rJl3^^^~2g^p — pcj 
/ _ 

= ±c v 7 h 2 c 2 + 2giP - f*c + 0 , 9 0 

V 2^//m 2 + 2gip 


+ 0(k 4 


0{k 4 


a 


H 

-io± 


H 




10 ±) 


(AS) 


a 


H 

-io± 


H 


= ±/?_ 


io±- 


Being different from a;± and uj± for a 0 and (3o, both o; G1 and o; G2 are gapful as long as 
gi > 0. In the nonrelativistic limit c —y oo, cu G1 diverges as well as oj±, and only the gapful 
mode with u; G2 survives. In the ultrarelativistic limit /x —> 0, a_i and (3-\ are independent 
of each other with the dispersion 

k 2 \ 


= ““ = ±C ('^ + 2VWp) + 0{k,) ' 


(A6) 


which remains gapful. Defining new operators 


F p i — 


0 0 0 

V 100 / 




Fp2 — 


0 0-j 
0 0 0 
i 0 0 


A 


(AT) 


/ 


we can write the modes corresponding to a_i and /3_i as e~ tFp2S i/j g = ^/p(l,0,s) r + 0(s 2 ) 
and e lFplS 'tfj s = y/p(1, 0, is) T + 0(s 2 ). Because we can obtain the nonmagnetic polar state 
V’p = y/PO-, 0, e l5 ) T j\[2 with F pl and F p2 as e -^p20osS-F pl smS)n/4,^ _ ^ we can re gard a_i 
and (3-i as the fluctuation from the ferromagnetic state to the polar state. 

In the case of g\ = 0, the number of NG modes changes as follows. In this case, two 
dispersions u; G1 and u; G2 become equivalent to uj± and u;± G for a 0 and /5 0 respectively, and 
a_i and (3 -1 also contribute to the type-II NG and Higgs modes. This is a consequence of 
the fact that the symmetry of the Lagrangian is enlarged from 17(1) x SO(3) to 17(3) and 
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broken generators for ip g included in u( 3) are F x> y iZiP i, P 2 , which have been just considered 
above. 

We finally refer to the linear-response theory, which gives the same results as those from 
the low-energy effective theory as well as other examples in the main part. The Bogoli- 
ubov equation can be obtained by substituting the fluctuations tp —» ip g + ue l ( k - x ~ ult + 5 6 7 ) -|_ 
v* e -i(k-x-ut+5) |-j ie or igj na ] Lagrangian flAll) . 


a; 2 3 4 /c 2 + 2pu — k 2 — F 

-G c u 


-G 


u 


2 /c 2 -2 pu-k 2 -F 


+ 0((u, v) 2 ) = 0 (A8) 


where F and G are given in the ferromagnetic ground state xp g as 


(<h ~ ,9i 0 0 ^ 


F = p 


^go - si o o'' 


, G = p 


V 


0 0 0 
0 0 0 


(A9) 


0 0 0 

y o 02 9l J 

We obtain the same dispersion relations as those discussed in the low-energy effective theory. 
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